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Structural and static properties of a classical two- 
dimensional (2D) system consisting of a finite number of 
charged particles which are laterally confined by a parabolic 
potential are investigated by Monte Carlo (MC) simulations 
and the Newton optimization technique. This system is the 
classical analog of the well-known quantum dot problem. The 
energies and configurations of the ground and all metastable 
states are obtained. In order to investigate the barriers and 
the transitions between the ground and all metastable states 
we first locate the saddle points between them, then by walk- 
ing downhill from the saddle point to the different minima, we 
find the path in configurational space from the ground state 
to the metastable states, from which the geometric proper- 
ties of the energy landscape are obtained. The sensitivity of 
the ground-state configuration on the functional form of the 
inter-particle interaction and on the confinement potential is 
also investigated. 

PACS numbers: 45. 05. -fx, 61.46. -fw, 73.22.-f 



I. INTRODUCTION 

Wigner suggested in 1934 that a liquid to solid phase 
transition should occur in a three-dimensional (3D) 
Fermi system at low densities 0. The quest for the 
observation of such a Wigner crystal has been the ob- 
ject of very intense and not stagnant work. After the 
first discovery of Wigner crystallization of electrons on 
the surface of liquid helium there has been con- 
siderable theoretical and experimental progress in the 
study of the mesoscopic system consisting of a finite 
number of charged particles, which are laterally confined 
by a parabolic potential and repel each other through a 
Coulomb potential. This system is the classical analog 
of the well-known quantum dot problem. These quan- 
tum dots are atomic-like structures which have interest- 
ing optical properties and may be of interest for single 
electron devices. These systems and their configurations 
have been observed experimentally, and are important in 
solid-state physics, plasma physics as well as in atomic 
physics. The classical approach is valid for quantum 
dots in high magnetic fields where the kinetic energy of 
the electrons is quenched, or for other classical systems, 
such as laser cooled ions in a trap |^ which are realized 
by electric and magnetic fields, trapped ions cooled by 
laser techniques ions in a radio- frequency (RF) trap 
(Paul trap) |^,^ or a Penning trap which can also 
serve as an illustration of 3D Coulomb clusters pO| , pl| . 
Very large Coulomb clusters have been created recently in 
strongly coupled RF dusty plasmas p2-p3 which are like 



a two-dimensional (2D) layered system. Examples of 2D 
Coulomb clusters are electrons on the surface of liquid he- 
lium fisll and electrons in quantum dots . The vortex 
clusters in an isotropic superfluid p^ , vortices in super- 
fluid He^ |8|jl9), vortices in a Bose-Einstein condensate 
stirred with a laser beam and in superconducting 
grains plf have many common features with those of 
2D charged particles Colloidal particles dissolved 

in water and placed between two glass plates are 

another example of an experimental system where clas- 
sical particles exhibit Wigner crystallization |2^]. Very 
recently, macroscopic 2D Wigner islands, consisting of 
charged metallic balls above a plane conductor were stud- 
ied and ground state, metastable states and saddle point 
configurations were found experimentally [ p6[ . 

In a finite system there is a competition between the 
bulk triangular lattice and the circular confinement po- 
tential which tries to force the particles into a ring like 
configuration. Those configurations were systematically 
investigated in Ref. |^ and a Mendeleev-type of table 
for these classical atomic-like structures was constructed. 
The spectral properties of the ground state configura- 
tions were presented in Refs. [ P8|j29| ] and generalized to 
screened Coulomb |30 3^ and logarithmic [ p0p^ -|3^ in- 
terparticle interactions. 

In the present paper we want to go one major step fur- 
ther and calculate not only all the different metastable 
states but also the saddle points between those local en- 
ergy minima and the path followed by the particles to 
transit between those energy minima. The present work 
is motivated by recent experimental work p6t where it 
was found that: i) some of the configurations did not 
agree with the previous theoretical published one, and 
ii) they were able to observe some of the saddle points 
which are the key configurations for transition between 
different stable (ground or metastable) states. Therefore, 
we also investigated the stability of the ground state con- 
figurations against the functional form of the confinement 
potential and the exact form of the inter-particle inter- 
action potential. 

The present paper is organized as follows. In Sec. 2, we 
describe the model system. In Sec. 3, our numerical tech- 
nique, used to obtain the ground and metastable states, is 
outlined. The technique we used to find the saddle point 
is similar to the Cerjan-Miller algorithm |^^. After the 
saddle points are found, we connect the saddle point to 
the global minimum or a local minimum by the 'walking 
downhill' method. Sec. 4 is devoted to the structural and 
static properties of the ground and metastable states for 

= 1 ~ 40. The configurations are analyzed and com- 
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pared with available experimental data and the results of 
previous theoretical approaches. The phase diagram for 9 
and 16 particles in the ground state with different func- 
tional forms of confinement potential and interparticle 
interaction is also calculated. The discussion on the sad- 
dle point is presented in Sec. 5, and the connecting path 
from the ground state to the metastable states is found, 
and we investigate the completely geometric properties 
of the energy landscape. Our conclusions are presented 
in Sec. 6. 



II. MODEL SYSTEM 

The model system consists of identical charged parti- 
cles interacting through a Coulomb repulsive interacting 
and moving in a 2D plane where they are confined by a 
parabolic potential 

The confinement potential Vijr') — ^frCuj^r^ is taken 
circular symmetric and parabolic, where m* is the ef- 
fective mass of the particles, q is the particle charge, 
wq is the radial confinement frequency and £ is the di- 
electric constant of the medium the particles are mov- 
ing in. Note that for the quantum dot problem an ad- 
ditional term appears in Eq. (1) which is the kinetic 
energy of the particles which is absent in our statical 
classical problem. Here the motion of the particles is 
restricted to the (x^y) plane. To exhibit the scaling 
of the system, we introduce the characteristic scales in 
the problem: rg = (2(7^/TOewQ)^/^ for the length and 
= {rnLo\(f l^e^Y^'^ foi' the energy. After the scaling 
transformations (r — > r/r^.E E/Eq), the Hamilto- 
nian can be rewritten in a simple dimensionless form as 

with = x'^ + and which only depends on the 

number of particles N. The numerical values for the pa- 
rameters wqj ^Oi E!q for some typical experimental systems 
were given in Ref. p7[| . 

III. NUMERICAL APPROACH 

Most of the previous works have treated the quantum 
mechanical problem of a small number of electrons. In 
the present paper, we consider only the classical system. 
Although a classical approach for the description of the 
behavior of electrons in quantum dots, in principle, is 
not applicable, it is possible that certain features of the 
classical system may survive in a quantum system. 



Due to the presence of confinement energy and 
electron-electron Coulombic interaction, a complete de- 
scription of the cluster system is complicated and can't 
be obtained analytically. The Monte Carlo simulation 
technique js^ is relatively simple and rapidly conver- 
gent and it provides a reliable estimation of the total 
energy of the system in cases when relatively small num- 
ber of Metropolis steps is sufficient. However, the ac- 
curacy of this method in calculating the explicit states 
is poor in certain cases. It becomes more difficult for 
clusters with a large number of particles, which have sig- 
nificantly more metastable states. To circumvent this 
problem we employ the numerical technique of Newton 
optimization which was outlined and compared with the 
standard Monte Carlo technique in Ref. ||28|1 . In this way, 
we are able to obtain not only the ground state but also 
metastable states. It also yields the eigenfrequencies and 
the eigenmodes of the ground state configuration. Now 
only a small number of calculation steps is needed to 
obtain the same accuracy. Moreover, using the modi- 
fied Newton approach, we can explore the stability of 
the system in its ground-state configuration through its 
spectrum. 

By studying the characteristics of the energy landscape 
and the energy barrier between the different local min- 
ima, we are able to find the saddle point configurations 
which are very important and are the key configurations 
for transition between different stable states. The tech- 
nique we used to find the saddle point is explained in 
more detail in Ref. [0, and is similar to the Cerjan- 
Miller algorithm |3^]. After the saddle points are found, 
we connect the saddle point to the global minimum or 
a local minimum by the 'walking downhill' method. In 
this algorithm the direction of the steepest gradient is 
followed to force the system to transit from the saddle 
point state to the local minimum state. Which minimum 
is finally reached depends on the initial step, therefore 
we repeat this procedure several times to determine both 
minima which the saddle point state connects. Thus the 
connecting path followed by the particles to transit be- 
tween those energy minima is found, from which the ge- 
ometric properties of the energy landscape are obtained. 

IV. GROUND STATE AND METASTABLE STATE 

In Table|we list for A^=l, 2,. ..,40 the energy per parti- 
cle E/N in the ground state and in the metastable states, 
where we also list the energy difference with the ground 
state AE/N. The configuration is indicated by the num- 
ber of particles in the different rings, the position of the 
center of the ring and the radius of the different rings, 
the width of the ring which is defined as the difference 
of the maximum radius and minimum radius in the same 
ring, and the energy of the lowest three normal mode 
frequencies of the ground state are also given in Table I. 
This table is rather exhaustive and should be compared 
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with a similar one published in Ref. for a logarithmic 
interacting system. 

For different values of N there exist different pos- 
sible values for E/N which are nothing else than the 
metastable states. The difference in energy between the 
metastable and the ground state is given in the third col- 
umn and the corresponding configuration in the fourth 
column. Note that with increasing N the number of 
metastable configurations increases and in general (but 
not always) the widths of the rings for metastable config- 
urations are larger and the central ring/particle is not ex- 
actly located in the center of the parabolic potential well. 
For sufficiently large iV, the simple ring structure gradu- 
ally disappears in the center and the triangular Wigner 
lattice appears. There is a competition between two 
types of ordering: ordering into a triangular-lattice struc- 
ture (Wigner lattice) and ordering into a shell structure, 
which leads to clusters with interesting self-organized 
patterns which show concentric shells at small N and 
hexagonal cores surrounded by circular outer shells at 
large N. 

The lowest non-zero normal mode frequency is a mea- 
sure for the stability of the ground state, it tells us how 
easy or difficult it is to deform this state. Therefore, in- 
tuitively we would expect that the value of this frequency 
would be correlated to AE/N, the energy difference be- 
tween the first metastable and the ground state. Those 
values are plotted in Fig. |^ as function of N. Notice that 
there exist such a correlation in general, but that this is 
not true for all TV- values, e.g. for N = 12, 18, 19, 20, 21, 30 
there is no correlation. 

The rings have sometimes a finite width which are 
shown in Fig. ^ as a function of N. Notice that the 
widths fall into three bands; i) width ^ 0.003 which is 
practically a perfect ring, ii) width ~ 0.02, and iii) width 
~ 0.5. Usually, but not always, the outer ring has the 
largest width. The width of the rings increases with in- 
creasing TV and at the same time the widest ring becomes 
often the next to outer ring. 

We compare our ground state configuration with avail- 
able experimental dat a [^6| and the results of previous 
theoretical approaches 30 3^j3^ . For very small num- 
ber of particles {N < 16), all theoretical and experimen- 
tal results for the ground state configurations are the 
same expect for iV=9 and 15 whatever kind of interpar- 
ticle interaction. The experimental observation for 
the ground state of 9 particles is (1,8) and for 15 parti- 
cles it is (4,11), which compares to our result (2,7) and 
(5,10), respectively. For 17 < < 30, the experimen- 
tal result and all the calculated patterns present three 
shells. Our result differs with the experimental data of 
Ref. Hi for N= 17, 20, 22, 24, 25, 27-30. Because of the 
discrepancy between some of the experimental configura- 
tions and the 'numerical exact' theoretical ground state 
configurations it is possible that experimentally the inter- 
particle interaction is not exactly a Coulombic potential 
and the confinement potential is not purely quadratic. 
Therefore, we investigated the effect of such deviations 



of these potentials on the ground state configuration. As 
an example we took A'^=9 and use confinement potentials 
y ~ r" and for the inter-particle interaction V ~ r~™. 
The resulting phase diagram is shown in Fig. ^. Notice 
that, depending on the values of n and m, the system 
can be either in the (1,8) or the (2,7) configuration. For 
the harmonic confined Coulomb interacting system i.e. 
(n, m) = (2,1) the system is in the (2,7) configuration 
but from the phase diagram it is clear that if we change 
the confinement potential slightly and make it more steep 
up to n ^2.2 the configuration (1,8) becomes the ground 
state. The experimentally determined ground state con- 
figuration for 9 particles was (1,8) [24]. 

There is also a difference with the experimental data 
and our results for 7V=16 particles. Therefore, we did 
the same investigation and present the phase diagram 
in Fig. H Notice that the harmonic confined Coulomb 
interacting system, i.e. (n, m) — (2,1) is again close to 
the phase boundary between the configuration (1,5,10) 
and (5,11). This is probably the explanation why the ex- 
perimental configuration pq] differs from our simulation 
results, since it is hard to guarantee that (n, m) is exactly 
(2,1) during the experiment. 

Notice that for both N=9 and N=16 the metastable 
configuration has an energy very close to the one of the 
ground state, the difference is less than 0.2%. These 
metastable configurations correspond indeed with the ex- 
perimentally observed ones. Consequently, an alterna- 
tive explanation for the difference with the experiment 
is that the experimental configuration got stuck in the 
metastable configuration. 



V. SADDLE POINTS 

Between metastable states and the ground state there 
are potential barriers. The system will prefer to trans- 
fer over the lowest potential barrier, which is the saddle 
point configuration between these energy minima, in or- 
der to transit from one stable configuration to the other. 
We plot in Fig. || the trajectories of the particles for the 
A^=5 system making a transition from the ground state 
(5) to the metastable state (1,4) and the saddle point 
connecting them. The trajectories of the particles can 
also be obtained by moving one of the particles to the 
center of the system. 

For 6 particles, the ground state (1,5) and the 
metastable state (6), corresponding to the hexagonal con- 
figuration, are obtained. Moreover, the unstable equilib- 
ria associated to saddle point configurations are also ob- 
tained, and the energy landscape is shown schematically 
in Fig. H There are two saddle points for this case, one 
of them is very close to the metastable state in both en- 
ergy and configuration, and will therefore be hard to see 
experimentally ||2^. In Fig. ||, the insets show the ar- 
rangement of the particles for the different states. Using 
the 'walking downhill method', we found the central par- 
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tide slowly moving to the periphery of the cluster. We 
would like to stress that the configuration with 6 parti- 
cles on a perfect ring is a saddle point state in contrast 
to the claim made in Ref. [Q. This can be understood 
from the following simple model calculation: if 3 particles 
are placed on a circle with radius A, on the corners of an 
equilateral triangle, and the other 3 particles on another 
equilateral triangle's corners with radius B rotated over 
60°, the energy is 

(3) 

where c = This function is shown in Fig. ^. It is 
clear that the perfect circle configuration i.e. c — ^ = 1 
is a saddle point, and that the minimum is obtained if 
3 particles move a bit to the center, and the other 3 
particles move away from the center (see the insets in Fig. 
0). Both shown metastable states are just connected by a 
rotation over 120*^. The two minima in Fig. ^ correspond 
to the same configuration in which inner and outer ring 
are interchanged. Comparing our results with the Fig. 
2 (W = 6: ground state, saddle point configuration and 
the hexagonal metastable state') of Ref. ||2^, we see that 
the other saddle point is observed experimentally. 

A list of the saddle point energies up to 20 particles is 
given in Table |l[ From this table, we notice that there 
is only one saddle point state for N = 3,4,5 particles. 
But, on the other hand it is well-known that there are 
(fc— 1) saddle points when there are k minima. For iV=3 
and 4 one saddle point is found, although there is no 
metastable configuration. The reason is that the saddle 
point state connects two equilateral ground state con- 
figurations which can be obtained from each other by a 
simple rotation. For the simple case of 3 particles, we 
show the energy surface and the corresponding configu- 
rations schematically in Fig. |[ Notice that there are 
always more saddle points than minima for N>6. With 
increasing the number of particles, more saddle point 
states are obtained and the energy landscape gets more 
complicated. For example for 9 particles, we obtain three 
saddle points and one metastable state. The results for 
the trajectories and energy landscape are shown in Fig. 
^. Again, the ground state configurations correspond- 
ing with the black and the white dot are connected by a 
simple rotation, i. e. a symmetry operation. 

VI. CONCLUSION 

We presented the results of a numerical calculation of 
the configurations of the ground and all metastable states 
and their energies, the system consisting of classical 2D 
charged particles that are confined in a parabolic con- 
finement potential for A^=l,...,40. These artificial atoms 



undergo configurational changes when the system tran- 
sits from the ground state to the different metastable 
states, or between the different metastable states. Such 
transitions move through the lowest energy barrier con- 
necting those states, i. e. through a saddle point. The 
connecting path from the ground state to all metastable 
states is found and the geometric properties of the energy 
landscape were discussed. 

Sensitivity of the configuration on the form of the con- 
finement potential and the interparticle interaction is in- 
vestigated and a phase diagram was obtained. This sen- 
sitivity on e.g. the form of the confinement potential is 
probably the explanation why the experimental configu- 
ration |2q| differs from our simulation results. 
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TABLE I. The ground state and the metastable states for 
A^=l,..., 40 Coulombic particles confined in a 2D parabolic 
well. We give the energies (E/N), AE/N, the shell structure 
{Ni,N2, ...), the radius and width of the shell, and the lowest 
three normal mode frequencies of the ground state configura- 
tion. 



TABLE II. The energies of the ground state, the 
metastable states and the saddle point states for different 
number of particles (A'^). 

Figure captions 

FIG. 1. The lowest eigenfrequency and AE/N as function 
of the number of particles. 



FIG. 2. The width of the different shells (logarithmic scale) 
as function of the number of particles. 



FIG. 3. The phase diagram for the ground state of 9 parti- 
cles. The dependence on the form of the confinement potential 
and the interparticle interaction is shown. 



FIG. 4. The phase diagram for the ground state of 16 parti- 
cles. The dependence on the form of the confinement potential 
and the interparticle interaction is shown. 



FIG. 5. The trajectories of the particles making a transi- 
tion from the ground state to the metastable state and the 
saddle point connecting them for 5 particles. 



FIG. 6. The energy landscape and transition between the 
ground state to the metastable states for 6 particles. 



FIG. 7. Part of the energy landscape and corresponding 
configurations near the metastable state for 6 particles. 



FIG. 8. Schematic view of the energy surface and projec- 
tion of the energy and the corresponding configurations for 3 
particles. 



FIG. 9. The energy landscape and transition from ground 
state to metastable states for 9 particles. 
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This figure "figl.png" is available in "png" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 06395 V 1 



M 
In 


P/M 


(Emetastable ~ 


L/OnTiguraDon 


riaulUS 
nf thp rinn 


\A/irl+h nf 
VVlUTn OT 

rinn 


Lowest Eigenfrequency 


1 


2 


3 


1 






1 












2 


0.75 


- 


2 


0.25000 





1.41421 


2.44949 




3 


1 .o lUo/ 




Q 
O 


u.4oD/y 


U 




1 . / o^uo 




4 


1 .oot)4o 




A 

4 




U 


1.25189 


1.41421 


1.86483 


5 


2.33845 


- 


5 


0.77948 





1 .02886 


1.41421 


1.94009 


2.36556 


0.0271 1 


1,4 


(0,0) 
0.98565 





6 


2.80456 




1,5 


(0,0) 
1.12182 





0.91889 


1.41421 


1.72935 


2.82476 


0.0202 


6 


0.94159 


0.13974 


7 


3.23897 




1,6 


(0,0) 
1 .25960 





1.09890 


1.41421 


1.45239 


8 


3.6689 


— 


1,7 


(0,0) 
1 .39768 





1.10688 


1.22689 


1.41421 


9 


4.08812 


— 


2,7 


0.18691 
1 .69865 


0.02355 
0.42157 


0.12681 


0.7754 


0.91669 


4.09426 


0.00614 


1,8 


(0,0) 
1 .53535 





10 


4.48494 


— 


2,8 


0.18189 
1 .82325 




0.47320 


0.0891 


0.97483 


0.9852 


4.48816 


0.00322 


3,7 


0.32295 
1 .99881 


0.06125 
0.39275 


11 


4.86467 


— 


3,8 


0.31407 
2.11186 


0.00303 
0.41407 


0.02451 


0.73727 


0.83817 


12 


5.23894 


— 


3,9 


0.31290 
2.22412 




0.42935 


0.53084 


0.89022 


1.15011 


5.24204 


0.0031 


4,8 


0.44386 
2.39909 




0.32645 


13 


5.60114 


— 


4,9 


0.43824 
2.50207 


0.00232 
0.36449 


6.002E-4 


0.75236 


0.75546 


14 


5.95898 


- 


4,10 


0.43499 
2.60686 


0.01827 
0.42444 


0.04940 


0.79446 


0.84578 


5.96269 


0.00371 


5,9 


0.56868 
2.77583 


0.01433 
0.2451 1 


15 


6.30758 




5,10 


0.56408 
2.87175 




0.30545 


0.45989 


0.68255 


0.75956 


6.31554 


0.00796 


1,5,9 


(-2.8 E-2, -6.7 E-3) 
0.83961 
3.04209 


0.07668 
0.41117 


16 


6.6499 




1,5,10 


(0,0) 
0.82994 
3.13165 




0.49686 


0.49237 


0.66241 


0.92179 


6.65235 


0.00245 


5,11 


0.55826 
2.97163 


0.00660 
0.43949 


17 


6.9829 


- 


1,6,10 


(0,0) 
0.94005 
3.39295 


0.01071 
0.27955 


0.05416 


0.54796 


0.55927 


6.98433 


0.00143 


1,5,11 


(-2.7 E-4, -5.3 E-3) 
0.81991 
3.22530 


0.01190 
0.58075 


18 


7.30814 


- 


1,6,11 


(1.9 E-3,-2.3 E-3) 
0.93020 
3.47888 


0.01259 
0.39138 


0.00614 


0.65492 


0.6782 


7.31522 


0.00708 


1,7,10 


(0,0) 
1.06145 
3.64611 


0.13829 
0.09771 


7.32316 


0.01502 


6,12 


0.68919 
3.31699 


0.42237 
0.61198 



1 



19 


7.63193 


- 


1,6,12 


(0,0) 
0.92663 
3.56465 




0.51668 


0.66759 


0.70531 


1.09892 




7.6328 


8.7E-4 


1,7,11 


(0.0001,0) 
1.04683 
3.72848 


0.00673 
0.21144 








20 


7.94961 


- 


1,7,12 


(0,0) 
1.03672 
3.81170 


0.01262 
0.28287 


1.031E4 


0.62728 


0.6926 




7.95623 


0.00662 


1,6,13 


(-0.0007, 0) 
0.91819 
3.65633 


0.01487 
0.82045 








21 


8.26588 


- 


1,7,13 


(0.0005, 0.0013) 
1.02997 
3.89626 


0.02028 
0.53913 


0.00317 


0.63738 


0.71452 




8.26645 


5.7E-4 


2,7,12 


0.15456 
1.27726 
4.05127 


0.02849 
0.44673 
0.61607 










8.26756 


0.00168 


1,8,12 


(0,0) 
1.16048 
4.04909 


0.40055 
0.44872 








22 


8.57418 


- 


2,8,12 


0.15569 
1.39275 
4.28533 




0.62765 
0.42346 


0.29341 


0.4018 


0.5569 




8.57568 


0.0015 


2,7,13 


0.15132 
1.26862 
4.13118 




0.43922 
0.74099 








23 


8.87758 


- 


2,8,13 


0.15121 
1.37540 
4.36583 


9.141 E-4 
0.55788 
0.60690 


0.12867 


0.4083 


0.58505 




8.87859 


0.00101 


3,8,12 


0.25859 
1.62401 
4.5251 1 


0.05289 
0.66983 
0.52010 








24 


9.17590 


- 


3,8,13 


0.25419 
1.60496 
4.60039 


0.00329 
0.48021 
0.62433 


0.02762 


0.40652 


0.5531 




9.17756 


0.00166 


3,9,12 


0.26352 
1.73670 
4.74997 




0.56748 
0.41325 








25 


9.47079 


- 


3,9,13 


0.25718 
1.71084 
4.82724 


0.01301 
0.51012 
0.48826 


0.11377 


0.5032 


0.52713 




9.47292 


0.00213 


3,8,14 


0.25286 
1.59306 
4.67414 


0.00985 
0.43153 
0.66009 










9.47485 


0.00406 


4,8,13 


0.36056 
1.83873 
4.83115 


0.02597 
0.42038 
0.61801 








26 


9.76273 


- 


3,9,14 


0.25454 
1.69487 
4.89949 


0.00247 
0.46781 
0.55831 


0.10409 


0.56681 


0.61503 




9.76383 


0.0011 


4,9,13 


0.35766 
1.93603 
5.05885 


0.04274 
0.42849 

0.25265 








27 


10.0509 


- 


4,9,14 


0.35397 
1.91746 
5.12748 


0.01222 
0.38138 
0.57653 


0.01311 


0.37880 


0.56307 




10.0527 


0.0018 


4,10,13 


0.35930 
2.04083 
5.27916 


0.07311 
0.50131 
0.46223 








28 


10.3356 




4,10,14 


0.35490 
2.01900 
5.34687 


0.03475 
0.46783 
0.33128 


0.05682 


0.17410 


0.47677 



2 





10.3378 


0.0022 


4,9,15 


0.35121 
1.90271 
5.19711 


0.01628 
0.35006 
0.58122 








29 


10.6181 


- 


4,10,15 


0.35181 
2.00193 
5.41431 


0.03407 
0.44070 
0.57886 


0.0391 1 


0.12706 


0.57988 




10.6193 


0.0012 


5,10,14 


0.46218 
2.23618 
5.57000 


0.07446 
0.43997 
0.41156 










10.6204 


0.0023 


4,11,14 


0.36293 
2.12559 
5.55935 


0.02363 
0.94989 
0.96329 








30 


10.8973 


- 


5,10,15 


0.45853 
2.21795 
5.63340 




0.32949 
0.48813 


0.33092 


0.35315 


0.43077 




10.8985 


0.0012 


4,11,15 


0.36734 
2.10775 
5.69516 


0.07206 
1.25661 
0.74043 










10.8999 


0.0026 


4,10,16 


0.35100 
1.98987 
5.48103 


0.05375 
0.45566 
0.67825 










10.9000 


0.0027 


1,5,10,14 


(-2.3 E-2, 1.7 E-2) 
0.68914 
2.45675 
5.78471 


0.07449 
0.62030 
0.58181 








31 


11.1739 




5,11,15 


0.45640 
2.31542 
5.84745 


0.04604 
0.39767 
0.45540 


0.02351 


0.11182 


0.27895 




11.1743 


4E-4 


1,5,10,15 


(0,0) 
0.68175 
2.43355 

5.84826 




0.51402 
0.67038 










11.1754 


0.0015 


1,5,11,14 


(-0.08,0.1) 
0.68601 
2.55743 
5.99410 


0.12595 
0.93520 
0.60146 










11.1756 


0.0017 


5,10,16 


0.45469 
2.20105 
5.69982 


0.01552 
0.35769 
0.57899 








32 


11.4466 




1,5,11,15 


(-1.3 E-2, 1.27 E-4) 
0.67693 
2.52730 
6.06078 


0.01588 
0.57582 
0.65479 


0.02971 


0.18804 


0.45106 




11.4479 


0.0013 


5,11,16 


0.45336 
2.29882 
5.90984 


0.02539 
0.42403 
0.56601 










11.4481 


0.0015 


1,6,10,15 


(-1.41 E-4, -1 E-2) 
0.77284 
2.64856 
6.06602 


0.03634 
0.41490 
0.61390 








33 


11.7156 




1,6,11,15 


(-3.2 E-3, -2.6 E-4) 
0.76976 
2.73747 
6.27609 


0.04088 
0.45135 
0.44249 


0.06805 


0.21724 


0.41934 




11.717 


0.0014 


1,5,11,16 


(4.6 E-3, 4.1 E-3) 
0.67207 
2.50792 
6.12124 


0.02433 
0.55453 
0.71356 








34 


11.9826 




1,6,12,15 


(0,0) 
0.77212 
2.83337 
6.47800 


0.03846 

0.60567 
0.48212 


0.23793 


0.35495 


0.78177 




11.9829 


3E^ 


1,6,11,16 


(5.4 E-3, -1.5 E-3) 











3 











0.76308 
2.71525 
6.33500 


0.01639 
0.42685 
0.60738 










11.9856 


0.003 


1,7,11,15 


(-4.4 E-3, -1.5 E-3) 
0.87023 
2.95104 
6.48557 


0.11285 
0.43662 
0.09455 








35 


12.2469 




1,6,12,16 


(0,0) 
0.76569 
2.80687 
6.53777 


0.00196 
0.54017 
0.53914 


0.06585 


0.46277 


0.47916 




12.2500 


0.0031 


1,7,11,16 


(2.2 E-3, -3.6 E-3) 
0.86132 
2.92300 
6.54591 


0.04942 
0.31919 
0.47110 










12.251 


0.0041 


1,6,13,15 


(-1.3E-2,-7.3E-3) 
0.77386 
2.93587 
6.67458 


0.18235 
1.14043 
0.64150 








36 


12.5108 




1,6,12,17 


(7E-4, 1.7 E-3) 
0.76095 
2.78682 
6.59542 


0.00271 
0.53950 
0.66210 


0.00895 


0.52837 


0.53908 




12.5111 


3E^ 


1,7,12,16 


(-2.9 E-3, 2.9 E-3) 
0.85806 
3.01116 
6.74952 


0.04137 
0.35782 
0.32542 










12.5124 


0.0016 


1,6,13,16 


(-2.5 E-3,-4.5 E-3) 
0.76430 
2.90429 
6.73797 


0.06647 
0.79964 
0.69310 










12.5135 


0.0027 


1,7,13,15 


(-0.041,-0.023) 
0.86813 
3.14480 
6.88005 


0.20913 
1.07946 
0.49131 










12.5139 


0.0031 


1,7,11,17 


(-2.31 E^, 3.8 E-3) 
0.85528 
2.90243 
6.60309 


0.03027 
0.29242 
0.71859 








37 


12.7719 




1,7,12,17 


(-3.3 E-5, 5.8 E-5) 
0.85232 
2.98912 
6.80497 


0.02693 
0.32128 
0.50272 


0.00321 


0.26342 


0.4995 




12.7724 


5E-4 


1,7,13,16 


(-2.21E-3,1.84E-4) 
0.85858 
3.10390 
6.94783 


0.06034 
0.59406 
0.43175 










12.7732 


0.0013 


1,6,13,17 


(0.013, -0.013) 
0.76718 
2.88302 
6.79140 


0.15906 
1.39644 
1.59430 








38 


13.0304 




1,7,13,17 


(-1.1 E-3, 1.6 E-3) 
0.85267 
3.07918 
7.00319 


0.02782 
0.49130 
0.52379 


0.00613 


0.16459 


0.46835 




13.0325 


0.0021 


2,8,12,16 


0.17423 
1 .22056 
3.39801 
7.13684 


0.22850 
1 .50692 
0.96666 
1.51773 










13.0327 


0.0023 


1,7,12,18 


(1 E-3, 4.1 E^) 
0.84689 
2.96973 
6.86200 


0.02664 
0.36496 
0.61396 









4 





13.0328 


0.0024 


2,7,12,17 


0.08029 
1.05842 
3.19102 
7.00656 


0.07477 
0.43210 
0.68153 
0.85628 










13.0331 


0.0027 


1,7,14,16 


(8.98 E^, 5.33E^) 
0.86305 
3.20222 
7.13834 


0.09294 
1.01212 
0.71661 








39 


13.2879 




2,8,12,17 


0.16869 
1.19687 

3.37497 
7.19592 


0.22552 
1.32722 

0.58608 
1 .79679 


0.22264 


0.31318 


0.36571 




13.2881 


2EA 


2,7,13,17 


0.13162 
1.06169 
3.28010 
7.20053 




0.48902 
0.97651 
0.67530 










13.2882 


3E-4 


1,7,13,18 


(9 E^, 0.0033) 
0.84868 
3.05915 

7.05762 


0.05232 
0.67493 
0.83314 










13.2885 


6E-4 


1,7,14,17 


(2 E-4, -5 E-4) 
0.85627 
3.17395 
7.19535 


0.04100 
0.86456 
0.82641 










13.2894 


0.0015 


2,7,14,16 


0.13300 
1.07631 
3.41100 
7.32553 


0.0028 
0.53723 
1.50092 
0.77203 








40 


13.5419 




2,8,13,17 


0.13522 
1.15191 
3.47914 
7.40261 


0.03091 
0.62836 
1.00580 
0.60815 


0.1242 


0.23757 


0.31601 




13.5423 


4E-4 


2,8,14,16 


0.14575 
1.16970 
3.61776 
7.51663 




0.74490 
1.44833 
0.69407 










13.5429 


0.001 


2,7,14,17 


0.12997 
1.06258 
3.37175 
7.39235 


0.01194 
0.46342 
1 .29363 
0.95848 










13.5434 


0.0015 


1,2,8,13,17 


(0.165, 0.214) 
0.39651 
1.27319 
3.47211 
7.39190 




0.61127 
1.56879 
1.70224 









E, r, CO are given in units of Eo, Tq, cOo/1/2 



Table I. Minghui Kong etal. 
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Number of 


Energy of 


Energy of 


Energy of 


Particles (N) 


Ground state 


Metastable state 


Saddle point 


1 


- 


- 


- 


2 


0.75000 


- 


- 


3 


1.31037 


- 


1.46201 


4 


1 .83545 


- 


1.92064 


5 


2.33845 


2.36556 


2.36829 


6 


2.80456 


2.82476 


2.82477 








2.82689 


7 


3.23897 


- 


3.27913 








3.28592 


8 


3.66890 


- 


3.68738 








3.68957 


9 


4.08812 


4.09426 


4.08813 








4.09530 








4.10628 


10 


4.48494 


4.48816 


4.48495 








4.48817 








4.48940 


11 


4.86467 


- 


4.87829 








4.87876 








4.88110 


12 


5.23894 


5.24204 


5.23955 








5.24209 








5.24368 


13 


5.60114 


_ 


5.61202 








5.61529 


14 


5.95898 


5.96269 


5.96335 








5.97235 


15 


6.30758 


6.31554 


6.30832 








6.31534 








6.31577 








6.31617 








6.31689 








6.31702 








6.31769 








6.32205 


16 


6.64990 


6.65235 


6.65117 








6.65490 








6.65547 








6.65620 








6.65648 








6.65706 


17 


6.98290 


6.98433 


6.98614 








6.98978 








6.98979 








6.99129 








6.99373 








6.99417 








6.99441 


18 


7.30814 


7.31522 


7.31522 






7.32316 


7.31535 








7.31700 








7.32046 








7.32095 









7.32154 
7.32390 


19 


7.63193 


7.63280 


7.63358 








7.63516 








7.64155 








7.64258 








7.64266 








7.64563 


20 


7.94961 


7.95623 


7.95637 
7.95638 
7.95639 
7.95640 
7.95641 
7.95654 
7.95693 
7.95695 
7.95709 
7.95725 
7.95739 
7.95817 
7.96316 



E is given in unit of Eo 
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